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I. Detection Theory and Detection Models 

In signal detection theory, the decision making portion 
of a detection system is called the and a detection 

experiment is the observation by a receiver of its input during a 
time interval. The input which is related to a target is called 
6-ignat, and the input which is not related to the target is called 
no^4e. In general, the observation is assumed to be of a known 
region which in some cases is called a Aziotut^on cztt. 

Signal detection theory is a basis for detection modeling. 

The detection models which are discussed here rely heavily upon 
it. When a detection experiment has been performed either the 

event = {At least one target was present in the region which 

was observed during the time of the observation.} or its comple- 
ment Hq will have occurred. In general, the models will specify 
that either the event D^^ = (The receiver decides at least one 

target was present in the region which was observed during the 

time of the observation.} or its complement Dq will have occurred. 
In terms of signal and noise, the events D^ and D^^ can be 
expressed as follows: Dq = {The receiver decides its input during 

the time of the observation was noise.} and D^^ = {The receiver 
decides its input during the time of the observation was signal 
and noise.} The theory which is the basis for detection models 
which have the above properties is called binary detection theory. 

Four events which are important in binary detection theory 
are indicated in the Venn diagram of Figure 1. 




Figure 1. Four events important to binary detection theory. 

The Venn diagram emphasizes a decision problem which is 
associated with a receiver. The problem is this; Under what con- 
ditions should the event occur, that is, under what conditions 

should the receiver decide that a target was present during the 
time of an observation? Detection theory may be able to provide a 
solution to this problem. If this is the case, the noise will 
generally be described as a random process and the signal will be 
described as either a deterministic or a random process. 

The following notation and terminology will be used: 

^ = P(D^|Hq) will be called the false alarm probability. 

^ = P(Dj^|Hj^) will be called the detection probability. P = P(Hj^) 
ull be called the prior probability. P is the probability that 
a. target will be in the region to be observed during the time of 
the observation. 

The input to a receiver is assumed to be a quantity whose 
square is proportional to power. The input at some time t^ will 
be symbolized by y(t^). The noise at t^ will be symbolized by 
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n(tj^) and the signal by s (t^) . The noise and signal processes 
will be represented by {n(t^),t^et} and {s (t^) , t^^et} where t 
represents the time interval during which the receiver 
observes the input. Often s(t^) and n(t^) can be assumed to 
be additive in which case Hq = {y(t^) =n(t^),t^et} and 
Hi = {y(t^) = n(t^) + s(t^) ,t^et}. 

Because of the finite quantity of information present at 
the input of a receiver, y needs to be measured at only a finite 
number of points in time in order to be adequately determined over 
an observation interval. For this reason, the noise and signal 
processes can be represented by {n(t^), i = l,...,m) and 
{s(tj^), i = l,...,m} where tj^,...,t^ are in the observation 
interval of length t. The noise process is then defined by a set 
of m random variables. To specify the noise process, one needs 
to specify only the joint distribution of the m random variables. 
If the signal process is deterministic, it is simply a finite set 
of values which is known before the detection experiment is per- 
formed . 
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II. 



Decision Criteria 



To simplify the discussion of decision criteria and decision 
rules, the stochastic process which represents a receiver's input 
will be assumed initially to be a single random variable Y. It 
can be called the decision random variable. The input process in 
this case is determined by the two conditional distribution func- 
tions FY(ylHQ) and FY(ylH^). 

The condition that a receiver's input is required to satisfy 
in order that will occur can be specified in terms of a deci- 

sion rule. For the assumed case, a decision rule is a rule which 
determines for every observable value of Y the decision that the 
receiver is to make. The decision rule can be considered to be a 
function <p(y) which relates the observable values y to the 
following two statements: 

d^: Decide the input was noise. 

dj^: Decide the input was signal and noise. 

Defining a decision rule cp (y) is equivalent to defining a set 
0 . such that D^ = {Y e fi} . 

The problem which was considered above can now be restated 
in the following way: What criterion should be adopted in order 

to determine a decision rule? A desirable characteristic for a 
criterion is suggested by the following argument: Consider the 

ratio of posterior probabilities 

P(H3^|y)/P(HQly) . 
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One might expect that the values of y in should make this 

ratio relatively large. But making this ratio large is equiva- 
lent to making the likelihood ratio L(y) large. This suggests 
that might be defined as follows: = (y: L(y) ^ K} where 

K is some positive constant which has yet to be specified. 

Four specific decision criteria are defined next in terms 
of n. For each criterion, ft has the above form with only the 
procedure for determining K being different. The decision cri- 
teria are: 

1. The Neyman-Pearson Criterion: Choose ft so that p^ is 

maximized subject to the constraint that p^ ^ a where a is 

a specified value. For a continuous decision random variable 
the constant K is chosen so that p^ = a. 

2. The Bayes Criterion: Choose ft so that the expected cost 

of a receiver's decision is a minimum. For a continuous decision 

random variable K = (^io~^OO^^Ol”^ll^ (1-P)/P if > Cqq and 

C-, > C, , where C. . is the cost of the event D. Pi H . . 

01 11 13 13 

3. The Ideal Observer Criterion: Choose ft so that the 

probability that the receiver makes an incorrect decision is a 
minimum. For a continuous decision random variable K = (1-P)/P. 

4. The Minimax Criterion: Choose ft so that the maximum 

expected cost of a receiver's decision is a minimum. For a con- 
tinuous decision random variable, K = (C^^q-Cqq/Cqj^-Cj^^) (1-P* ) /P* 
if C^Q > Cqq and C^^^ > Here P* is the value of P which 

would make the expected cost of a receiver's decision a maximum 

if a Bayes decision rule cp„ were used. 
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A more general discussion of the above criteria would involve 
the notion of a randomizing rule. 

If a model is adopted which specifies the conditional dis- 
tributions and a decision rule, then the values of p^ and p^ 
can in principle be determined. The pair of values (p^,p^) is 
called a receiver operating point. If the decision rule results 
from one of the four likelihood ratio criterion listed above, then 
it will involve the parameter K through the relation 

n = {y: L(y) ^K}. And, for a given value of K, since determines 

(Pf,Pd)f a unique operating point result. By varying K, a set 

of operating points can be generated which is called an ROC curve. 

Different ROC curves can be produced by changing one or both of 

the conditional distributions. Note, a change in a conditional 

distribution function implies a change in the signal or the noise. 

Decision rules which result from applying the four likelihood 

ratio criteria in a model in which the input process is determined 

by a set of m random variables can be expressed in terms of a set 

as follows: fi={(yT,...,y): L(y,,...,y)^K} where K is 

1 m 1 m 

specified in the same way as it is in the corresponding case in 
which the input process is determined by a single random variable. 
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III. Two Detection Models 



A particular detection model will now be examined. In 
the model, the receiver's input process is defined by; 



Hq; {y (t^) = n (tj^) , i = l,...,m} and 

{y (t^) = n (t^) + s (t^) , i = 1, ,m} 

where the values of y are measured every At units of time. 

The noise process {n(t^), i = l,...,m} is assumed to consist 
of a set of m independent normal random variables each with 
mean zero and variance . The signal process {s(t^), i = l,...,m} 
is assumed to be deterministic. Thus, the input process consists 
of m independent normal random variables Yj^,...,Y^ each with 
variance and each with mean zero when the target is not 

present and with mean s^ = s(t^) when the target is present. Note, 
such a model might be used to obtain an optimistic estimate of a 
detection systems performance, since all the information about 
the signal is assumed to be known. 

The result of the application of a likelihood ratio decision 
rule in the model can be expressed in terms of a random variable 
Z. This random variable is called the cross correlation statistic 
and it is defined by Z = Y s.Y.. However, it is more convenient 

-i ‘-11 

to express the results in terms of a random variable V which is 
defined by V = Z/a^. In terms of this random variable. 



and 



p^ = 1 - $(v*) 
p^ = 1 - f(v*-dJs) 



7 



where $ symbolizes the standard normal distribution function and 
V* = (1/0,) (o^ In K+ (1/2) Is.), dh = I s?/o,- 

^ -L J. 

The parameter d is called the detection index and, since 




Recall, the input stochastic process is assumed to repre- 
sent a quantity whose square is proportional to the input power. 
Therefore, the average power input to the receiver over the time 
interval t during which the receiver observes the input from a 
region can be approximated by I y?/m. The average signal power 
can be approximated by S = I s?/m and the average noise power 
by I n?/m. And the expected average noise power can be approxi- 
mated by N = I a^/m = . in these terms, d = m(S/N) . 

The receiver's bandwidth will be represented by BW. 

Assuming m corresponds to the value determined by the sampling 
theorem, that is, m = t/At = 2t(BW), the detection index can be 
written as d = 2t(BW)(S/N) or, defining Nq = N/BW/ as d = 2t(S/NQ) 
A model in which the signal process is not deterministic 
will be considered next. The model is defined by 

Hq: {y^^ i = 1/ ... /in} 

{y^ = nj^ + Sj^, i = l,...,m}. 

Here the noise process is a set of m independent normal random 

variables with mean zero and variance as before. However, 

now the signal process is also a set of m independent normal 

random variables with mean zero but with variance a^. The result 

s 

of applying a likelihood ratio decision rule in this model can be 



8 



expressed in terms of a random variable X which is defined by 

X = ^ Y?. When a target is not present, the statistic X/a^ has 

a chi-square distribution with m degrees of freedom and when a 

target is present, the statistic X/(a^+a^) has a chi-square 

distribution with m degrees of freedom. Hence, p^ = 

and p, = P [x^^x*/ (a^+a^ ) } where x^ is a chi-square random 
d m s m 

variable with m degrees of freedom and x* is a number which 
is determined by the decision rule. 

When the target is not present, since X/a^ has a chi-square 
distribution, the variance of X is 2mo‘' and the mean is ma^. 
When the target is present, since X/(a^+a|) has a chi-square 
distribution, the variance of X is 2m(a^+a|) ^ and the mean is 
m(a^+a|). By the central limit theorem, as the number of degrees 
of freedom m of a chi-square random variable becomes large it 
can be approximated by a normal random variable with the same mean 
and the same variance. Hence, for sufficiently large m, p^ 
and p^ can be approximated as follows: 

p^ = 1 - $ [ (x*-ma^ ) / ( 2ma‘* ) ^] 

and 

p^ ^ 1 - ${ [ (x*-ma^ ) - ma|] / (2ma‘* ) ^] 

where (a^+a|)^ has been approximated by o'*. This approximation 
can also be written as: 

p^ ~ 1 - $ (v*) 
p^ = 1 - $(v*-d*^) 

with V* = (x*-ma^)/ (2ma‘*) ^ and d = (ma|) ^/2ma‘* . 
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Note, over the observation interval t, ^ s?/m 
the average signal power so S = 

the expected average signal power and, as before 
approximately the expected average noise power, 
model the detection index can be written as 

d = j (S/N) ^ = t(BW)(S/N)2. 



is approximately 
approximately 
N = is 

Hence, for this 
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IV. Detection Model Applications 

When a likelihood ratio decision rule was used with either 
of the two models discussed above, the following result was obtained 

Pj = 1 - $(v*) 

h 

p^ = 1 - $ (v*-d ) 

In the case of a known signal (deterministic process) , a cross 
correlation receiver is required, and d = 2t(BW)(S/N). This will 
be called Case I. In the case of a Gaussian signal, a square law 
receiver is required and, for S/N << 1 and (BW) large, 
d=t(BW) (S/N)^. This will be called Case II. In both Case I 
and Case II, the detection index d is a function of (S/N). 

The two models can be related to radar and sonar systems, for 
example, by using the radar equation or sonar equation to deter- 
mine (S/N) . 

In some radar and sonar detection models, it is assumed 
that the decision rule used results in a required p^ and in 
addition that a minimum value of p^ is specified in the sense 
that a target is said to be detectable only when p^ is greater 
than or equal to the minimum value. This minimum value of p^ 
along with the required p^ define what could be called a mini- 
mum acceptable signal-to-noise ratio sonar models, this 

minimum acceptable signal-to-noise ratio determines the detection 
threshold DT by the relation DT = 10 log (S/N) If the mini- 
mum acceptable value of p^ is .5, then DT is usually called 
the recognition differential RD. 
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A passive sonar detection model will now be discussed. It 
an example of a model which uses the concept of a minimum value 
p^. It also illustrates a method of dealing with non-station- 
r noise and signal processes. In the model, the event detection 
the event {X^g ^ 0} where Xg^ represents signal excess, a 
ndom variable with expected value SE. The source level, noise 
vel , directivity index, and the detection threshold or recogni- 
on differential are also random variables with expected values 
! NL, DI , and DT respectively. The passive sonar equation 
n be written in terms of these random variables or in terms of 
?ir expected values. In the latter terms, SE = FOM - TL, where 
M SL - (NL-DI) - DT is the figure of merit. To determine the 
obability of a detection P(X„„^0), the distribution of the 

bhi 

idom variable X„„ must be specified. 

This model could be considered to be equivalent to one in 
' ch the event {Xg^^O} is the event that p^ is greater than 
linimum acceptable value. Although this definition may appear 
re reasonable than the one which defines {X„„^0} as the event 

bili 

tection, it would likely be no more useful for developing tac- 
s, designing exercises, or testing detection devices. 

The noise and signal processes should be effectively sta- 
onary during a "look" at a resolution cell, that is, over an 
-egration period if a Case I or Case II model is used to relate 
p^, and the signal-to-noise ratio. If either of these models 
used, the signal-to-noise ratio for an integration period will 
the value of a random variable, and, therefore, p^ will be 
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conditioned on the random variable signal-to-noise ratio. The 
receiver is usually assumed to be readjusted so that remains 

relatively constant from "look" to "look" at a given resolution 
cell in applications of the model. 
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V. Search Theory and Search Models 

Search theory deals with problems such as that of determining 
the optimal allocation of search effort in a region or that of 
determining the probability of detecting a target within a given 
length of time for a particular search plan. In treating these 
problems, false alarms are often ignored. In effect, the time 
spent on false alarms is assumed to be small relative to the time 
scale of the problem, however, the cost associated with a false 
alarm is not assumed to be negligible. 

In general, it is assumed here that during a search a 
target will move through resolution cells surveyed by the detec- 
tion system of a searcher. In the model of such a search, the 
probability that the input from a resolution cell containing the 
target will cause the receiver to declare the target is present 
is p^. In systems with operators, the operator can delay the 
decision regarding the presence or absence of a target in a reso- 
lution cell and can recall the input of adjacent resolution cells. 
Such detection systems are difficult to model. To illustrate 
this, consider the following hypothetical system consisting of an 
active sonar system plus an operator. Assume that p^ as a func- 
tion of the signal-to-noise ratio has been determined for the 
operator in a laboratory experiment by forcing the operator to 
respond after a single "look" at a resolution cell. Operationally, 
the probability that the operator detects a target will not be 
given by the laboratory determined value of p^ if the operator 
delays the detection decision for several "looks." Various models 
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have been proposed to deal with this situation. In one, a "three- 
out-five" detection criterion has been adopted. The criterion 
states that an operator will declare a target if, out of five 
consecutive inputs from adjacent resolution cells, at least three 
are such that they would have caused the operator to declare a 
target present in a single "look" experiment. An input which 
satisfies this condition is referred to as a success. With this 
criterion, the probability that a target will be detected with 
m consectuive "looks" at the target is found as follows: Deter- 

mine the 2 ’'''' sequences of successes and failures of length m 
which could result with m consecutive "looks." The probability 
is equal to the sum of the probilities of occurrence of those 
sequences for which the three-out-of-f ive criterion is satisfied. 
The probability of a particular sequence will depend on the value 
of p^ and 1 each of the m "looks" in the sequence. 

Unfortunately, in most situations, the number of sequences to be 
considered is too large to be computationally tractable even if 
the probabilities in the sequence can be determined. 
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VI. The Probability of Detection During a Search 

The problem of determining the probability that a target 
will be detected in a search at or before n "looks" have occurred 
is generally basic to the solution of search problems. Let N 
represent the ntimber of the "look" at which detection first occurs, 
then this probability can be expressed as 

n 

P(N^n) = I P(N=i). 
i=l 



It can also be expressed as 



n 

p(N^n) = 1 - ”TT (1-g^) 

i=l ^ 

where g^^ = P(N=i|N^i-l) is the probability that the event detec- 
tion occurs on the i— look conditioned on the event detection 
has not occurred earlier. The second expression for P(N^n) is 
generally of greater interest than the first one, since g^ can 
be directly related to operational parameters such as target range. 

■f~ H 

Note, if the resolution cell being examined on the i— look is 
empty, then g^ will be zero. 

A continuous analog to P(N^n) can be developed as follows: 

Consider a model in which the time for a look is At, and if 

detection occurs on the i— "look," the time of detection is 

th 

iAt. The probability of detection at or before the n — "look" then 
can be expressed in terms of time by noting that P (N^n) = P(T^nAt) 
where T represents the time of detection. Now define the detec- 
tion rate to be y (nAt) = (1/At) P [T=nAt 1 T^ (n-1) At] . The probability 
P(T^nAt) can be approximated in terms of a related quantity y(t) by 
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P(T^t) = 1 - exp{- j Y(T)dt} 

0 

where Y(t) is related to Y (nAt) by the requirement that 
P(T^t) equals P(N^n) for t = nAt. In some models where it 
would be useful to use P(T^t), the continuous approximation to 
P(N^n)/ it might not be reasonable to model the detection process 
in terms of single "look" probabilities. If this is the case, 

Y (t) must be defined by some other means. The visual detection 
model developed in the OEG Report No. 56 is an example of this. 

Generally, in search operations, the target will be moving 
on a relative track C. In some cases, it may be useful to con- 
sider the track to consist of segments ,€2 / . . • . Let 
represent the event the target is first detected on track segment 
. The probability of first detecting a target while it moves 
along C can then be expressed as follows: 

n 

P(T) = 1 - exp{- I F(C. ) } 

i=l ^ 

where F(C^) is called the sighting potential on the track seg- 
ment and exp{-F(C^)} = P(r^|fj^_ 3 ^n H . It is conven- 

ient to use the concept of sighting potential in the analysis of 
some search types, for example, "ladder searches". 

The track of a target which moves in a plane can be described 
by equations x = x(t) and y = y(t). These equations describe 
the relative track of the target with respect to the searcher 
(often called the observer) if x and y refer to a coordinate 
system in which the searcher is located at the origin. In this 
case, r the range of the target is given by r = [x^+y^] 
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If in a model the detection capability of a searcher 
against a target is assumed to be independent of the time at which 
the target occupies a particular position on the target's track, 
then time can be replaced by target distance along the track. 



f w 
•'0 



(T)dT 



This can be done by using the transformation equation s 
where w is the target's speed relative to the searcher and the 
zero of time has been chosen to coincide with the start of the 
search. If, in addition, the detection capability of the searcher 
against the target is assumed to depend only on the range of the 
target, then the sighting potential for a target's track C for 



a continuous model can be expressed in terras of either t or s 

ft rS 

Y [r (T) ]dx = 

0 •' 0 



as follows: F(C) = 

dt/ds = l/w(s) . 



Y [r (a) ] da/w (a) since 
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VII. Experimental Validation of Detection Models 

Consider a detection system which could be satisfactorily 
modeled by a continuous detection rate function that depends only 
on r, the target's range, if a suitable detection rate function 
could be found. Suppose on some basis that a particular detec- 
tion rate function is proposed? An estimate of the suitability 
of the proposed function could be made using standard statistical 
tests if appropriate operational data could be obtained. For 
example, the number of detections at various ranges for straight 
line encounters. If such data were available, its agreement with 
the cumulative distribution function of the target's range at 
detection deduced from the proposed detection rate function could 
be determined by use of the Kolmogorov-Smirnov test. This technique 
could be applied even though the number of encounters in which a 
target is not detected was unknown. In this case, the ciamulative 
distribution function to be compared is Fj^(r|D) where D = 

{detection} and R represents the target's range at detection. 

To derive F (rlD) or equivalently the density function 
f„(r|D) given a continuous detection rate function, one can pro- 

r\ 

ceed as follows: First find the joint density function which de- 

scribes the distribution of the target's rectangular coordinates 
(X,Y) at its detection point. To do this, set f^^ |D)AxAy = 

fy (x,y)AxAy/p where p = P(D). This can be done, since 

A / X 

[ (Y=y) Od] = (Y=y) . That is, Y can be considered to take on a 

value only when a detection occurs. Note, p is the value of the 

integral of f^ ^(x,y) over all pairs of values (x,y) . Now note that 
X # Y 
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^(x,y) = f„ I „ (y I x) (x) and that the conditional density 

X / Y Y I X X 

f I (y|x) can be written as 
Y I X 



fyix^yl^^ = (1/w) y[r (y) ]exp{-j Y [r (n) ] dn/w} . 

— 00 

Given (x) and Y(r), (r | D) can be found by first trans- 
forming to r and a by using the transformation equations x = r 
cos a and y = r sin a and then finding the marginal distribution 
of R as follows; 



fR(r) 



2-n 

(r , 0 ) da , 

0 



An example of the above procedure will now be given. Suppose 



f^(x) = 



J 



l/2a 



X ^ a 



X > a 



and y(r) = 



kh/r^ 

0 



X ^ a 



X > a 



The transformation to r and a gives 

r sin a 

I (l/ 2 awp) (kh/r^)exp{- 

I 



f^^^(r,a|D) = 



kh (r ^cos^+n^ ) ^"^^dTi/w} 
I r cos a I ^ a 



I r cos a I > a 



For large a, this can be simplified by letting a -> «>, after 

f ^ 

noting that p = l/2a p(x)dx = W/2a where W is the sweep 
width. In the limit, ^ f„ „(r,a|D) = (kh/wWr ^ ) exp {-kh/wr^ 1-sin a)}, 
and after performing the integration over a, 

fj^(r|D) = (W/2r2) {l-$(kh/wr^)^}. 

Note, which emphasizes the fact that in 

general detection need not take place during a straight line encounter. 
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VIII. Sweep Width Determination for a Random Track Angle 

Now consider the problem of determining probability p of 
detecting a target during an encounter which could be described 
by a straight line encounter model if the angle cp between the 
target's track and searcher's track were known. Suppose it can 
be assumed that lateral range and the relative speed 
w = [u^+v^-2uv cos cp] ^ are values of random variables whose joint 
distribution is known. The probability p of detecting a target 
during a straight line encounter can then be expressed as 



P 



fCO 



— oo 



I p(x,w) f^^ ^(x,w)dxdw 
0 



where p(x,w) = p (det |x = x 0 = w) and f^^ ^(x,w) is the joint 
distribution of the lateral range X and the relative speed 
If X and 'F are independent random variables and if 
p(x,w) as a function of w is adequately described by a linear 
approximation, then p can be approximated as follows: 



where 

w = E(i|)) 



P 



^OO 

p (x,w) fj^(x)dx 

— oo 



wf ^^(w) dw 

h 



•2 IT 

(u^+v^-2uv cos cp) f^((p)dcp. 

0 



An average sweep width can be defined by 



W = 



p(x,w)dx. 



If $ is uniformly distributed so f^(<p) = w is the value 

of an elliptic integral of the second kind. An average sweep width 
based on this value of w should be used in the Random Search 
Formula when either u or v are not appropriate. 
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IX. A Parallel Sweep Search Model 

A model of a parallel sweep or "ladder" search will be 
considered next. In the model, the target's velocity is zero 
and its position is uniformly distributed in a rectangular region. 
The searcher moves along a series of n parallel tracks to that 
n straight line encounters will occur in a complete search of 
the region. The track spacing is s and the tracks are numbered 
from 1 through n. The search geometry is shown in Figure 2. 



n 



I 

ns 



Figure 2. Search geometry for a parallel sweep model 



-> z 



If the target's horizontal coordinate is a random variable 
Z in the coordinate system shown in the figure, then P(det|z=z) 



n 



= l-exp[- 5! F(Xj)] where Xj = z - (j-%)s is the target's 



j=l 



• th 



lateral range from the j — track and 

/•ns 



P(det) = 



P (det Z=z) dz/ns . 
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In some cases, it is convenient to use the following coordinate 
system: Let I be the random variable which is the number of 

the track closest to the target. Relabel the tracks so that the 

number of the track closest to the target is 0, then = -js + x 

where now j = -K, . . . , 0 , . . . ,M where n = M + K + 1 and the value 

of K and M are determined by the value of I, that is, by 
the location of the target. Now 



m 

P(det| X = X n I=i) = l-exp[- I F(x.)] = p(x,i) 

j=-k ^ 

and 



Hence, 



P ( de 1 1 1 = i ) 



fS/2 

p (x,i)dx/s = p(i) . 

-s/2 



n n 

P(det) = I p(i)P(I=i) = (1/n) I p(i) 
i=l i=l 



or 



n [s/2 



P(det) = (1/ns) I 

i=l J 



p(x, i) dx 



-s/2 



If n is large and the detection law is such that y (r) is 

essentially zero for targets at ranges equal to or greater than 

one or two track spacings away, then for almost all values of i 

and j, p(i) will be approximately equal to p(j) and it can 

fS/2 

be assumed that p(det) = (N/Ns) p(x;s)dx where 

-s/2 

oo 

p(x;s) = l-exp[- I F(-js+x)]. 

j=-co 
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To illustrate an application of the above result, suppose 
Y (r) = kh/r^. Then F(Xj) = 2kh/wXj . If it is assumed that the 
number of tracks and the track spacing is such that p(i) = p(j) 
for almost all values of i and j , then 
m “ 

I Fix.) - I F(-js+x) = (2kh/w) (Tr^/s^)csc^ (TTx/s) 
j=-k ^ j=0 

so 

p(x;s) = 1 - exp{- (2khiT^/ws^ ) csc^ (irx/s) } . 

and 

P(det) = 2 ${(tt/2)^ (W/s)} - 1 

h 

where W = 2(2irkh/v) is the searcher's sweep width against the 

target. The ratio W/s is called the coverage factor. A plot of 
P(det) is given as a function of W/s in the National Search 
and Rescue Manual. 

The track length which is expended in an unsuccessful search 
is nb by the assumptions of the parallel sweep model. Since 
the area searched is nsb, the probability of detecting a target 
with a random search in the area with track length nb is given 
by P(det) = 1 - exp{-W/s}. This value could be compared to that 
given by the parallel sweep model in order to obtain a relative 
measure of effectiveness for a parallel sweep search. 
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X . The Optimal Allocation of Search Effort Problem 

The problem of determining an optimal allocation of effort 
for a search is discussed next. Here, an allocation of effort is 
an optimal allocation if it maximizes the probability of detecting 
the target. 

The class of searches which will be considered are those 

which can be described by a random search model that satisfies the 

following conditions: The target is fixed in a region of area A. 

The region consists of m subregions which are determined by the 

th 

conditions that in the i— subregion the sweep width is a con- 
stant w^ and the probability density function is a constant 

where A^^ is the area of the subregion and p^ is 
the target's prior probability of being there. (The subregions 
are assumed to be niimbered so that ^ ^ ^ 

target is detected in a subregion, the searcher must be in the 
subregion when it is detected. The probability of detecting the 
target is then 

n 

P(det) = I {l-exp(-(p- ) }p-A. 
i=l 

where (p . = VI. i./ A. is called the search effort density. 

^1 X 1 1 

Now consider the problem of determining the solution set 

(0 . which maximizes P(det) subject to the constraints 

1 op 
^ n 

Z = I (A./W.)<p. and (p . ^ 0 . Note, I is the available track 
i=l ^ ^ 

length. This is a nonlinear programming problem for which the 
solution set is: 
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i in 



InW. p. -d/X) I (A./W^)ln p. + Z/X 
^ ^ 3 J D 



<P 



1 op 



not in ^ 



where X = ^ ^ “ {1,2,. ..,k} where k is determined 

by the condition that if k + 1 were included in fi then 0, 

Suppose the model described above were used to determine 

an optimal allocation of track length for a search. The result 

would be a set of values Z> = (A./W.'jw. which maximize 

1 op i' 1 op 

the probability of detecting the target. Given detection will 
occur, in what order should the subregions be searched in order 
to minimize the expected time until the target is detected? The 
following procedure will effectively minimize the expected track 
length to detection and gives a suitable order assuming the 
searcher's speed remains constant so that minimizing the expected 
track length to detection is equivalent to minimizing the expected 
time to detection. 

Divide Z into minimum increments AZ consistent with 
the random search model and then allocate the increments in the 
following order: Allocate AZ to the Region 1. This will be 

assumed to be effectively consistent with the prescription deter- 
mined by the formula for cp . given AZ. If this search is 

unsuccessful, determine the optimal allocation prescribed for 
2AZ. This will effectively indicate whether the next increment 
AZ should be assigned to Region 1 or Region 2. Continue in this 
fashion until detection or until Z has been expended. Note, 
this is an approximate procedure in the sense of the model. 
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For a discussion of the above procedure, let E(l|l^) 1) 
symbolize the expected track length to detection given detection 
will occur by the time track length I has been used. That the 
above search procedure will effectively minimize E(l|ls£) for 
an optimal allocation of track length is suggested by the follow- 
ing argument: Divide I into segments Ail as indicated above, 

then the probability that the target is detected on or before the 
t h 

i — step of the search P(L^iAJl) is, in effect, a maximum, 

since no other allocation of search effort iA£ will give a 

significantly greater value for any value of i. Hence, 

{ iAX, |L:^il) = P (L^iAS. I is also effectively a maximum for any 

value of i since P (L^iAil | L^£) = P (L^iAil) /P (L^sil) and P(Ls;Jl) 

must be equal to its maximum value for any alternate procedure. 

n 

Therefore, E(l|l^£) = ^ [ 1-F (iA£ | L^£) ] is effectively a 

i=l ^ 

minimum for the procedure. 

The above procedure is the one which would have been followed 
if after each unsuccessful search one calculated new prior proba- 
bilities using Bayes procedure and then did the search prescribed 
for a track length A£ with the new priors. For a discussion of 
this point, see OEG Report No. 56. An application, in effect, is 
discussed in Reference 4. 
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XI. Target Position Probability Distributions Which Change in Time 
Consider the following model in which a target moves with 
constant course and speed in a plane with its position specified 
in a rectangular coordinate system: The target's course and speed 

are independent random variables $ and U respectively 
whose density functions and given. The joint 

density function f^^ ^{x,y;t) of the target's coordinates is 
needed, but only f „(x,y;0) the joint density function for time 

A , X 

zero is given. 

For such cases^ found in principle as 

A , I 

follows: First it is convenient to use polar coordinates and with 

X = p cos a and y = p sin a and to set 

f{p,a;t) = f^^Y( p cos a, p sin a;t) . 

Note, to first order, Pftarget is in AA at t) = f(p,a;t)AA. 




Figure 3. Problem geometry. 
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Referring to Figure 3 



/•«> r27T 



f (p,a; t) AA = 



f [r{u,(j)) ,6(u,(i)) ;0]AA'f^((|))fy(u)d4)du 



0 0 



where the functions r(u,(J)) and 9(u,(})) are such that if the 
target's polar coordinates were r(u,<J>) and 6{u,<})) at time 0, 
then with the speed and course u and <J) they would be 
(p,a) at t. The functions can be expressed implicitly by 

r^ = p^ + (ut)^ - 2p ut cos (a-({)) and r cos 9 = p cos a - ut cos cj> . 

The integration is, in effect, a sum over all possible pairs of 
values of u and (j). Each pair determines a vector ut which 
translates target positions at time 0 to target positions at 
time t such that AA is AA' translated without rotation or 
distortion. So, AA' = AA, independent of u and (}) . Therefore, 



f (P/Ot; t) 



f [r(u,<})) ,9 (u,(J)) ;0]f^((J))fy{u)d(f)du 



0 0 



As an example, suppose 





with f*(4)) = , 0 < (b ^ 2 tt, but with U = u known. In this 

$ 2ir 



case , 




r^ (u,(b) 
2o^ 



and, using r^ = p^ + (ut) ^ - 2p ut cos (a-<))) 
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f(p,a;t) = 



2iTa ' 



2n - [ p^+ (ut) ^-2p ut cos (a-(f) ) ] 



d4> 

2tt 



2ira 



1 2a 

T e 



[p^+(ut)2] 



2tt 



f2iT cos (a-((>) 

e ° d(l> 



J 

0 



hr [p^+(ut) 



2tto^ 



2a'^ ..r- • j ^ ciput'i 

r^J 



'^of- 



where Jq is the zero order Bessel function of the first kind 
and i = To obtain f ^(x,y;t) , substitute P^ = 

A , X 

in the above result. In this case, since y(x,y;t) is symmet- 

ric about the origin, the marginal distribution of the target's 
range R from the origin is of interest. In general. 



~ cos a,p sin a;t) = pf(p,a,-t) 



and 



fj^(P;t) = 



2tt 



fR,e (P,a;t)da 



so, in this case. 



[p' + (ut)^] 



f^(p;t) =^e ^ 



'0^ a- 



The distribution function F (p;t) for this case is given in 

K 



Reference 5. 
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